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Abstract 

Interior  ballistic  performance  of  unslotted^- 
slngle-perforated  stick  propellants  in  a  large- 
caliber  guq  could  be  influenced  significantly  by 
the  eroslveppurning  effect.  To  analyze  erosive^ 
burning  phenomena,  a  comprehensive  model}  with 
special  emphasis  on  the  Interaction  of  turbulence 
and  combustion} has  been  formulated.  Experimental¬ 
ly,  a  test  rig  using  a  center-perforated  cylindri¬ 
cal  propellant  grain  with  large  web  thickness 
(rl.Ocm)  was  designed,  constructed , and  tested. 
The  location  of  the  instantaneous  propellant 
burning  surface  was  determined  by  the  use  of  realn 
time  X-ray  radiography.  Based  upon  X-ray  Images, 
the  Instantaneous  burning  rate  of  NOSOL-363  stick 
propellant  was  found  to  be  much  higher  (up  to  3.2 
times)  than  the  strand  burning  rate  under  strong 
cross-flow  conditions.  / e  ,  ,  4-.j  "  \ 


Nomenclature 


a 

A1 

b 


C,..Cn 


CP 

Cr,.Cr2 

cs 

cu 

DR1 


DR2 

Ea 

Eb 

F 

H 

Ir 

Jr 


Flux  model  absorption  coefficient,  a-1 
preexponential  factor  of  1th  species, 
m3/kmol  -s 

covolume  of  Noble-ADel  equation  of 
state,  m3/kg 

turbulence  constants  for  k  and  e 
equations 

constant  pressure  specific  heat,  J/kg-K 
turbulence  constants  for  transport 
equation  of  reactedness 
specific  heat  of  stick  propellant,  J/kg-K 
turbulence  constants  for  eddy  viscosity 
group  of  species  pyrolyzed  from 
propellant  surface  having  delayed 
reactions 

delayed  reaction  species  generated  from  0 
and  F  species 
activation  energy,  J/kmol 
black-body  emissive  power,  •  oT4,  J/wt-a 
fuel-rich  species  pyrolyzed  from 
propellant 

"IYih<  *  uiui/2,  stagnation  enthalpy  of 
1th  species,  J/kg 

outward  radiation  flux  In  positive 
radial  direction,  W/m2 

Inward  radiation  flux  in  negative 
radial  direction,  W/m2 
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turbulence  kinetic  energy,  m2/s2 
interaction  length  of  X-ray  and 
propellant  sample 

linear  attenuation  coefficient  for 

photoelectric  absorption 

oxidizer  species  pyrolyzed  from 

propellant  surface 

pressure,  N/m2 

final  product  species 

radiative  heat  flux  absorbed  by  the 

solid  propellant  surface,  W/m2 

surface  heat  release  due  to  pyrolysis  at 

reference  temperature,  J/kg 

propellant  burning  rate,  m/s 

inner  radius  of  perforation,  m 

outer  radius  of  stick  propellant,  m 

gas  constant,  J/kg-K 

reactednesses  defined  in  Eqs.(3A)  and 
(3B) 

relectlvllty 

flux-modeling  scattering  coefficient,  1/m 
time,  s 

temperature,  K 

gas  velocity  in  axial  direction,  m/s 
gas  velocity  in  radial  direction,  m/s 
molecular  weight  of  1th  species, 
kg/kmol 

axial  coordinate,  m 

mass  fraction  or  1th  species,  i  could 
represent  F,  0,  DR1 ,  DR2,  or  P 


Greek  Symbols 

t  turbulence  dissipation  rate,  m2/s3 

ca  surface  emlsslvlty  of  solid  propellant 

*  thermal  conductivity,  W/m-k 

u  dynamic  viscosity  of  gas,  N-s/m2 

ut  turbulent  viscosity,  N-s/m2 

vi  number  of  kmoles  of  ith  species 

p  density,  kg/m3 

a  Stefan -Boltzmann  constant,  W/n^-K4 

u>i  rate  of  production  of  1th  species, 

kg/m3-s 


Subscripts 

4  centerline  value 

eff  effective 

g  gas-phase 

1  internal  perforation  region 

p  propellant 

s  surface 


Olaorltlcala 


mass-weighted  average  quantity 
Superscripts 

"  fluctuation  quantity  in  Favre  averaging 


Introduction 

The  Interior  ballistic  performance  of 
unslotted  single-perforated  stick  propellants  in  a 
large-caliber  gun  could  be  influenced 
significantly  by  the  effect  of  erosive  burning. 
In  order  to  analyze  and  predict  the  performance  of 
large-caliber  gun  systems  using  stick  propellants, 
the  erosive-burning  characteristic  of  the 
propellants  must  be  determined.  In  the  past, 
erosive  burning  of  various  solid  propellants  was 
studied  extensively  under  different  cross-flow 
situations.1  However,  prior  to  this  study,  no 
suitable  theoretical  model  was  specifically 
developed  to  simulate  erosive-burning  processes 
occurring  inside  the  stick  perforation  of  low 
vulnerability  ammunition  (LOVA)  propellants.  This 
is  due  to  the  fact  that  LOVA  propellants  have  a 
distended  flame  structure  and  introduce  stronger 
interaction  of  combustion  and  turbulence.2  This 
interaction  cannot  be  properly  modeled  by  previous 
theoretical  approaches.  Experimentally, 
instantaneous  burning  surface  locations  and 
erosive-burning  rates  along  the  internal 
perforation  of  propellant  grain  could  not  be 
determined  because  of  the  limitation  of 
measurement  techniques  and  hostile  experimental 
conditions.  It  is  apparent  that  a  more 
sophisticated  theoretical  model  and  more  advanced 
measurement  technique  is  needed  to  characterize 
the  erosive-burning  processes  occurring  inside 
perforations  of  stick  propellants.  The  objectives 
of  this  study  are:  1)  to  formulate  a  theoretical 
model  to  describe  erosive-burning  phenomena,  with 
special  emphasis  on  turbulence-combustion 
interaction;  2)  to  demonstrate  the  feasibility  of 
using  real-time  X-ray  radlography3  to  measure  the 
instantaneous  burning  surface  locations  of  an 
axisymmetrlc  stick  propellant  grain;  3)  to  display 
images  obtained  from  erosive-burning  tests  showing 
the  instantaneous  web  thickness  of  stick 
propellants;  and  4)  to  discuss  advantages  and 
limitations  of  the  real-time  X-ray  method  for 
erosive  burning  measurements. 

Theoretical  Modeling 


The  physical  model  considered  in  this  study 
consists  of  two  different  regions:  the  gas  phase 
and  the  solid  phase.  The  gas-phase  region  is 
occupied  by  an  axisymmetrlc  turbulent  reacting 
flow  inside  the  perforation  of  an  unslotted  stick 
propellant  (see  Fig.  1).  The  solid  phase  is  a 
NOSOL-363  propellant  with  an  initial  geometry  of  a 
thick-walled  cylindrical  tube  with  uniform  inner 
and  outer  diameters.  To  study  the  erosive  effect, 
only  the  internal  perforation  surface  is  allowed 
to  burn. 


shorter  than  the  characteristic  times  associated 
with  conduction  in  the  solid  propellant  and 
surface  regression.  This  assumption  has  usually 
been  adopted  in  past  erosive-burning  studies. 
Second,  it  is  assumed  that  the  heat  flux  in  the 
direction  normal  to  the  propellant  surface  is 
dominant  because  of  strong  heat  feedback  from  the 
flame  zone.  Thus,  the  axial  heat  conduction  is 
considered  to  be  negligible. 

As  pointed  out  by  Flfer,2  the  combustion 
mechanism  of  nitramine-based  LOVA  propellants  is 
very  similar  to  that  of  nitrate  ester  homogeneous 
propellants,  such  as  nitrocellulose/nitroglycerin 
(NC/NC)  double-based  propellants.  Following  the 
description  of  flame  structure  by  Wu  et  al.2  and 
Kuo  et  al.,3  the  reaction  regions  for  gas/solid 
interface  and  gas  phase  consist  of  five  different 
zones,  i.e.,  preheat,  foam,  fizz,  dark,  and 
luminous  flame  zones. 

Under  the  conditions  of  zero  or  low 
cross-flow  velocities,  Kubota°  found  significant 
flame  stand-off  distances  for  both  double-base  and 
nitramine  propellants.  Under  strong  cross-flow 
conditions,  the  flame  could  be  highly  distorted  or 
even  form  pieces  of  flamelets;  hence  the  location 
of  the  heat  release  zone  produced  by  the  reaction 
of  pyrolyzed  ingredients  could  be  quite  different 
from  the  location  of  pyrolysis.  Finite-rate 
chemical  reaction  kinetics  and  the  interaction  of 
turbulence  and  combustion  must  therefore  be 
considered  in  the  erosive-burning  model  in  order 
to  simulate  combustion  processes  inside  the 
perforation  of  a  stick  propellant  made  of  nitrate 
ester  homogeneous  propellants. 

In  modeling  chemical  reactions,  the  following 
reaction  mechanism  was  adopted  from  Wu  et  al .  *• 
For  surface  pyrolysis  of  stick  propellant,  we  have 

Stick  propellent  *  vg0  •  *rF  •  vM)D*1  (  ;  ) 

where  0  represents  the  oxidizer-rich  gases  (such 
as  NOg)  pyrolyzed  from  the  surface,  F  represents 
the  fuel-rich  species  (such  as  CHgO) ,  and  the 
remaining  species  are  designated  as  the  first 
group  of  delayed  reaction  species,  DR1  .  Three 
chemical  reactions  considered  to  take  place  in  the 
gas  phase  are 


V  *  vo°  * 


(2A) 


(2B) 


Fig.  1  Flow  and  temperature  fields  inside  an 
unslotted  single-perforated  stick  propellant  grain 


W**  *  V'  (2C) 


where  DR2  represents  the  intermediate  product  due 
to  the  reaction  in  the  fizz  zone.  After  a  delay 
in  the  dark  zone,  the  second  group  of  delayed 
reaction  species  reacts  in  the  luminous  flame  zone 
to  generate  the  final  product.  Croup  1  of  delayed 
reaction  species  also  forms  the  final  product  ir, 
the  final  flame. 


There  are  two  basic  assumptions  in  the 
formulation.  First,  the  mean  flow  is  treated  as 
quasi -steady  since  the  flow  residence  time  la  much 


2 


The  stoichiometric  coefficients,  vj'g>  in 
Eqs.  (2A-2C)  are  determined  from  the  mass  balance 


of  the  above  three  reactions.  The  heat  of 
formation  of  these  representative  groups  of 
species  is  obtained  from  the  chemical  equilibrium 
code  (CEC765.7  Detailed  procedures  are  described 
in  Refs.  4  and  5. 

It  is  well-known  that  one  major  aspect  in  the 
modeling  of  turbulent  reacting  flow  lies  in  the 
handling  of  the  time-averaged  chemical  production 
term  which  is  complicated  by  the  interaction  of 
turbulence  and  combustion.  The  effect  of 
Inevitable  fluctuations  of  flow  properties  in  this 
source  term  is  difficult  to  model.  One  logical 
way  to  circumvent  this  difficulty  Is  to  adopt  a 
two-variable  Joint  probability  density  function 
(pdf)  to  obtain  statistically  averaged  physical 
properties.® 

The  approach  used  in  modeling  the  chemical 
production  term  is  similar  to  that  of  Janlcka  and 
Kollmann^  in  their  studies  of  turbulent  Hg-alr 
diffusion  flames.  The  present  approach  differs 
in  the  selection  of  Independent  variables. 
Janlcka  and  Kollmann  used  the  mixture  fraction 
with  a  single  reactedness;  the  present  approach 
adopts  two  reactednesses.  This  is  due  to  the  fact 
that  in  premixed  flames  it  is  impossible  to  use 
mixture  fraction  as  one  of  the  Independent 
variables.  One  specific  feature  of  the  present 
approach  is  similar  to  that  of  Bray  and  Moss,10 
who  also  used  a  Joint  pdf  with  two  reactednesses 
as  Independent  variables.  In  their  studies,  two 
consecutive  reactions  were  used  to  simulate 
hydrocarbon-air  combustion,  while  a  set  of 
competitive  and  consecutive  reactions  are 
considered  in  the  present  study. 

The  two  Independent  variables  chosen  for  the 
pdf  are  reactednesses  of  reactions  (2A)  and  (2B) . 
They  are  defined  as 


<3A> 


,  ,  TPti  *  *0*1. u 

*0*1 .b  '  rD*l.u 


<3B) 


The  transport  equations  for  mean  reactednesses  and 
mean  square  fluctuations  of  reactednesses  are 
given  as 


applied  to  turbulent  diffusion  flame  by  Janlcka  et 
al.114  This  approach  was  not  adopted  in  the 
present  study  because  of  the  extraor d 1 nar y 
numerical  effort  required  to  solve  the  evolution 
equation.  The  second  approach,  to  assume  a 
functional  relationship  between  the  pdf  and 
selected  variables,  has  been  used  extensively  in 
calculation  of  different  turbulent  diffusion  or 
premixed  flames. 9. 10  In  this  study,  a  simple  beta 
functional  form  was  chosen  to  relate  the  pdf  to 
the  two  reactednesses  and  their  variants.  The 
beta  pdf1®  is  defined  as 


where 


f>  <a>  - 


fi  -  X)8'1 
l(i  -  M8'1** 


(6) 


i] 


(7) 
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Under  the  assumption  that  R *  and  Rg  are  not 
correlated,  the  Joint  probability  function,  ,®j(Ra. 
Rg) ,  can  be  expressed  as 


#>j(*A.  «,)  (9) 


By  using  pdf,  mean  chemical  production  terms  in 
Eqs.  (4)  and  (5)  can  be  expressed  as 

^  -[  f  w*»  (10) 

J  0  J0 


*/0I0  <*  -  (11) 


where  the  Instantaneous  chemical  production  term 
is  given  by 


Jx  >r  r  »r  5c  »i  * 
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.Vxp(-E,  £  b  .  Yru,  (12a) 
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where  R  can  be  either  R^  or  Rg. 

In  general,  there  are  two  approaches  in 
determining  the  relationship  between  the  pdf  and 
the  two  selected  variables.  One  approach  is  to 
solve  the  equation  of  evolution  of  the  pdf.  The 
modeling  of  the  pdf  evolution  equation  was 
proposed  by  Doparo,11  Pope,12  and  O’Brien, '3  end 


FO*  *  ’  S 

For  premixed  flames,  in  circumstances  in 
which  all  species  are  assumed  to  have  the  same 
coefficient  of  diffusion  transport,  the  carbon 
element  mass  fraction,  Zc,  and  the 
Schvab-Zeldovlch  variable,  £ ,  can  be  regarded  as 
constants,  that  is. 


3 


Zc  -  - 


Turbulent  constants  used  in  the  above  governing 
equations  are  listed  in  Table  I. 


e  -  t0  -  (-2-fiuF  -  con»t*nt 


equal  to  one,  i.e., 


Table  I  Constants  used  in  turbulence  modeling 


all  species  ls 

Constant 

Cl 

c3  c3  CA  cu  CRV 

CR2 

(15) 

Value 

1.0 

1.3  1.57  2.0  0.09  2.7 

1.79 

The  mass  fractions  of  all  five  species  can  be 
solved  from  five  algebraic  equations  [including 
Eqs.(3A),  <3B),  (13),  (I'D.  and  05)]  after  Hx 
and  are  solved  from  Eqs.(if)  and  (5). 

The  turbulent  flow  Inside  the  center 
perforation  of  the  stick  propellant  can  be 
described  by  the  following  set  of  conservation 
equations  for  mass,  momentum,  total  enthalpy, 
transport  equation  of  turbulent  kinetic  energy, 
and  transport  equation  of  turbulent  dissipation 
rate,  as  well  as  the  equation  of  state  in  the 
Favre-averaged  form.  The  following  equations  were 
derived  with  the  following  assumptions:  1)  Favre- 
averaged  flow  properties  are  quasi-steady;  2) 
boundary  layer  is  axisymmetrlc;  3)  body  force  and 
gas-phase  radiation  heat  transfer  are  negligible; 
V)  Soret  or  Dufour  effects  are  small;  and  5) 
Pick's  law  of  mass  diffusion  is  valid. 


jj  (t*  0  »)  -  JJ-  (r*  0  . 


1  “  It  • r  *  &  •  ii  {?  <**».«  &  •  Jf 


To  determine  the  Instantaneous  temperature 
profile  in  the  stick  propellant,  as  well  as  the 
surface  temperature  distribution,  a  transient 
one-dimensional  heat  conduction  equation  is 
considered.  This  equation  is  given  as 


•[<vv +  <w] 


which  takes  into  account  the  subsurface  radiation 
absorption  for  translucent  propellants  such  as 
NOSOL-363.  The  determination  of  the  source  terms 
containing  both  outward  and  inward  radiant  fluxes 
requires  solution  of  the  flux -transport  equations. 
Based  upon  a  two-flux  model,  the  two  transport 
equations  for  Ir  and  Jr  can  be  given  as 

4(rl„)  1 

-  -(*+») rlr  ♦  arE^  +  Jf  ♦  f  *r(Ir+Jr>  (25) 
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A  set  of  boundary  and  initial  conditions  must 
be  specified  to  complete  the  theoretical 
formulation.  At  the  solid-gas  interface,  balances 
of  mass  and  energy  fluxes  can  be  written  as 
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where  m  -  1  for  axisymmetrlc  flews. 

Turbulent  viscosity,  uf  is  expressed  in 
terms  of  k  and  e  as 


In  deriving  the  above  equations,  the 
following  relationships  were  used. 
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where  QSlpef  ls  defined  as  the  net  surface  heat 
release  at  reference  termperature ,  Ts  refi 
(qpad)  etls  the  net  radiation  flux  on  the 
interPIce,  and  ls  related  to  Ir  and  JP  by  the 
following  equation. 
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Other  boundary  conditions  at  the  solid-gas 
interface  are 


(30) 
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3)  treatment  of  subsurface  radiation 
absorption  with  two-flux  equations. 


which  represent  non-slip  condition,  mass 
continuity,  and  zero  fluctuation  of  reactednesses , 
respectively.  The  boundary  conditions  for  the  k 
and  c  equations  at  the  near-wall  region  are 
Identical  to  those  of  Arora  et  al.'5 

When  the  flow  Inside  the  perforation  of  the 
stick  propellant  Is  not  fully  developed,  flow 
properties  are  considered  to  be  either  continuous 
or  constant  at  the  edge  of  the  boundary  layer, 
i  .e. , 


f-r.  is  -  is  .  o 
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When  the  flow  is  fully  developed,  the  symmetric 
conditions  are  used  at  the  axis,  l.e., 

(33) 
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Experimental  Approach 


Experimental Apparatus 


Parallel  to  the  theoretical  study,  an 
experimental  investigation  has  been  carried  out  to 
observe  the  erosi ve-burning  phenomena  and  to 
establish  a  data  base  for  model  validation.  Since 
stick  propellants  used  in  large-caliber  guns 
generally  have  cylindrical  geometry  with  single 
perforation,  it  is  advantageous  to  use  a  test 
sample  with  cylindrical  geometry  to  attain  close 
simulation.  Under  normal  circumstances,  phenomena 
occurring  in  a  cylindrical  grain  cannot  be 
observed  by  conventional  flow  visualization 
techniques.  This  is  why  most  researchers  in  the 
past  have  used  two-dimensional  propellant  slabs  as 
test  samples.  Becuase  a  real-time  X-ray  is  now 
available  for  erosive-burning  measurements,  a 
teat  propellant  grain  with  cylindrical  geometry 
can  be  adopted.  In  order  to  allow  a  relatively 
long  test  duration  (a  few  seconds)  for  surface 
regression-rate  measurement,  a  thick-walled 
cylindrical  grain  with  web  thickness  in  the  order 
of  1  cm  has  been  used.  A  set  of  NOSOL-363 
propellant  samples  with  2.54  cm  outer  diameter  and 
0.467  cm  inner  diameter  has  been  processed  for 
erosive-burning  tests. 


The  initial  condition  for  the  solid-phase  heat 
.onductlon  equation  is 


Tp  ‘  V  »T 


Boundary  conditions  for  Eq.(24)  and  the  radiative 
heat-flux  equation  at  the  gas-solid  Interface  are 
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and  at  the  outer  radius  are 
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The  distinct  features  of  the  model  presented 
:n  this  paper  are: 

1)  consideration  of  gas-phase  reactions  with 
pdf  closure  to  account  for  the 
interaction  of  combustion  and  turbulence; 

2)  use  of  two  reactedness  parameters  to 
simulate  delayed  reactions  in  the  dark 
zone  of  homogeneous  propellants;  and 


A  test  rig  was  designed  and  constructed  to 
provide  two  different  configurations,  A  and  B  as 
shown  in  Figs.  2  and  3,  respectively.  The  test 
chamber  of  configuration  A  is  made  of  stainless 
steel  with  thermocouple  ports  for  subsurface 
temperature  measurements  at  various  locations . 
The  main  body  of  the  test  chamber  of  configuration 
B  is  made  of  LEXAN  materials  to  facilitate  X-ray 
penetration.  Configuration  B  is  also  equipped 
with  a  fast-actuating  ejection  mechanism  using  an 
explosive  bolt  to  achieve  Interrupted  burning  and 
propellant  recovery.  Ignition  of  the  propellant 
charge  in  the  driving  motor  is  accomplished  by 
activating  the  electric  primer  via  a  current  input. 
The  combustion  products  of  the  propellant  charge 
flow  from  the  driving  motor  into  the  test  section 
through  a  converging  graphite  nozzle.  As  hot 
gases  enter  the  stick-propellant  perforation,  the 
test  sample  is  heated  quickly  to  reach  a  fully 
ignited  condition.  The  combustion  product  gases 
generated  from  the  propellant  sample  Join  the 
product  gases  from  the  driving  motor  and  flow 
through  the  exit  nozzle  to  the  ambient.  Ir. 
-■onf iguratlon  A,  the  propellant  grain  can  also  be 
recovered  by  replacing  the  exit  nozzle  with  a 
bursting  diaphragm  assembly. 


■SCTTrta  •  it  A  JftXVISKvm 


Fag.  2  Schematic  diagram  of  test  rig  assembly  for 
studying  erosi ve-burning  phenomena  of  stick 
propellants  (Configuration  A) 
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Fig.  3  Schematic  diagram  of  test  rig  assembly  for 
studying  erosive-burning  phenomena  of  stick 
propellants  (Configuration  B) 

In  configuration  B,  some  physical  processes 
occur  Inside  the  test  rig.  Instead  of  measuring 
subsurface  temperatures ,  the  surface  regression 
processes  are  observed  and  recorded  using 
real-time  X-ray  radiography  techniques.  Employing 
these  two  configurations,  the  erosive-burning  test 
rig  can  provide  the  following  special  features: 

1)  the  instantaneous  inner  surface  location 
of  the  test  sample  can  be  observed  by 
means  of  a  X-ray  radiography  system; 

2)  the  driving  motor  can  generate  high- 
pressure  high-velocity  product  gases  to 
flow  over  the  propellant  sample  surface; 

3)  burning  can  be  confined  to  the  Internal 
surface  by  using  flame-retardant  grease 
on  the  external  surface; 

A)  ejection  of  the  exit  nozzle  and 
propellant  grain  Into  a  water  tank  can  be 
achieved  by  Incorporating  an  explosive 
bolt  assembly  with  a  timing-deal y  device; 
and 

5)  the  test  chamber  pressure  can  be  altered 
by  using  interchangeable  exit  nozzles. 

Data  Acquisition  System 

During  the  test  firing  with  configuration  B, 
the  instantaneous  burning  surface  location  of  the 
propellant  grain  is  filmed  by  a  real-time  X-ray 
radiography  system.  Figure  *1  shows  the  layout  of 
various  components  in  the  radiography  system.  A 
continuous  X-ray  is  generated  from  the  Phillips  MG 
321  constant  potential  X-ray  system.  Two  X-ray 
tubes  [one  with  a  focal  spot  combination  of  0.2  x 
0.2mm/3.0  x  3.0  ram  (MCN  167-160  kV) ,  the  other 
with  a  focal  spot  combination  of  1.2  x  1.2  mm/A.O 
x  A.o  mm  (MCN  321*320  kV)]  are  included  in  the 
system  to  satisfy  different  penetration  depth  and 
spatial  resolution  requirements.  A  lead  diaphragm 
is  installed  at  the  exit  port  of  the  X-ray  tube 
head  to  limit  the  angle  of  divergence  of  the  X-ray 
beam  and  to  confine  the  beam  to  the  measuring 
section  of  the  erosive-burning  test  rig  (this  also 
reduces  unnecessary  radiation  exposure  to  the 
surroundings).  A  second  lead  diaphragm  with  a 
larger  opening  is  placed  Jn  front  of  the  image 
lntenslfler  to  reduce  scattered  X-ray  radiation 
and  resulting  noise  on  the  fluorescent  screen. 


Fig.  A  Layout  of  real-time  X-ray  radiography 
system 


After  passing  through  the  test  rig.  X-ray 
signals  are  transformed  to  fluorescent  light 
signals  on  the  output  screen  of  a  tri-field  image 
intensifler  (Precise  Optics,  Model  PI  2A00  ATF , 
A",  6",  or  9"  field  diameter).  The  input  fluor  of 
the  image  intensifler  is  made  of  cesium  iodide 
with  a  decay  time  constant  of  650  ns,  and  the 
output  fluor  is  a  p20  type  with  a  85  ns  decay-time 
constant.  These  time  constants  are  short  enough 
to  allow  the  motion  analysis  system  to  operate  at 
its  maximum  framing  rate  without  generating  image 
blur. 

The  fluorescent  light  signal  output  from  the 
image  intensifler  is  recorded  by  a  Spin  Physics 
2000  Motion  Analysis  System.  This  system  can 
record  up  to  2000  fully  digitized  frames  per 
second,  or  up  to  12,000  digitized  pictures  per 
second  with  adjustable  playback  speed.  Two 
cameras  can  be  used  simultaneously  to  film 
different  views  of  the  event.  The  motion  analysis 
system  consists  or  the  following  subsystems: 

a)  a  camera,  with  solid-state  image  sensor. 
The  picture  information  goes  to  the 
console  from  the  camera,  and  is  processed 
into  a  frequency-modulated  carrier  that 
is  recorded  on  tape. 

b)  a  main  electronic  bin,  which  contains  the 
record  and  playback  electronics,  along 
with  video  output  circuitry. 

c)  a  tape  transport,  which  drives  the 
half-inch  tape  cassette  at  a  maximum 
speed  of  250  Inches  per  second. 

d)  a  Dage  camera  and  Panasonic  VCR,  which 
record  the  event  at  30  frames  per  second. 

Digitized  data  is  stored  on  a  high-intensity 
magnetic  recording  tape  and  can  be  transferred  t' 
the  digital  image-processing  system  frame  by  frari 
for  analysis  through  an  IEEE-A88  interface.  The 
digital  image-processing  system  consists  of 
several  major  components: 

a)  a  Quantex  (QX-9210)  digital  image 
processor  with  two  pipeline  point 
processors,  each  with  a  random  access 
image  memory  of  A80  x  6A0  x  12  bits.  The 
processors  perform  real-time  image 
enhancements,  Including  noise  reduction, 
image  subtraction,  arbitrary  contrast 
control,  roam,  and’  zoom.  The  processor 
can  also  be  used  to  conduct  such  high¬ 
speed  image  analyses  as  brightness 
histogram,  local  contrast  stretch,  area 
and  point  brightness  measurements, 
calibrated  length  and  area  measurements, 
Sobel  edge  enhancement,  and  so  forth; 

b)  a  mass  storage  device,  consisting  of  both 
a  2. A  MB  dual  8"  floppy  diskette  drive 
and  a  50  MB  hard  disk  for  the  storage  of 
image  data;  and 
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c)  a  high-resolution  videograph  printer  for 
producing  large  high-quality  pictures  on 
hard  copies. 

The  advantages  of  using  real-time  X--ay 
radiography  for  erosive-burning  measurements  are 
summarized  below: 

1)  No  viewing  window  is  needed  to  record  the 
erosive-burning  event  since  it  is 
possible  to  observe  the  surface 
regression  process  occurring  in  a  steel 
or  optically  opaque  chamber. 

2)  X-ray  radiogrphy  is  a  non- lnt rus i ve 
tehcnique  which  will  not  affect  the 
erosive -burning  processes. 

3)  The  technique  can  be  applied  to  either 
2-D  or  axisymmetrlc  propell ant-gr al n 
geometry. 

it)  It  is  a  real-time  data  recording  system. 
Pictures  captured  in  the  test  event  can 
be  played  back  immediately  after  the  test 
firing  at  a  lower  speed  for  visualization 
and  analysis. 

5)  All  data  are  in  digital  form  and  are 
convenient  for  recording,  analysis, 
transfer,  and  storage. 

Data  Reduction 

Despite  the  advantages  listed  above,  the  data 
reduction  procedures  for  deducing  the 
Instantaneous  internal  diameter  of  a  burning 
propellant  sample  is  not  straightforward. 

Theoretically  speaking,  an  ideal  radiography 
image  of  a  hollow  cylindrical  solid  propellant 
with  known  geometry  can  be  determined  from  the 
assumptions  that  1)  X-ray  radiation  is  generate- 
from  a  point  source,  via  an  infinitely  small  focal 
spot,  2)  X-ray  is  only  attenuated  by  photoelectric 
absorption,  and  3)  distribution  of  X-ray  Intensity 
over  the  input  screen  of  an  image  lntenslfler  is 
uniform.  As  shown  in  Fig.  5.  the  X-ray  Intensity 
distribution  over  the  input  screen  of  an  image 
lntenslfler  can  be  evaluated  by  the  following 
“quatlon: 

1(9)  -  Ig  expC-m  L( 8) ]  (39) 

IOEAL  X-RAY  POINT  SOURCE 
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CYLINDRICAL 
PROPELLANT  SAMPLE 


LO. — J  ' — INPUT  SCREEN  OF 

IMAGE  INTENSIFIES 

Fig.  5  X-ray  intensity  distribution  of  a  hollow 
cylinder  under  idealized  conditions 


wh-rre  m  is  the  linear  attenuation  coefficient  :  or 
photoelectric  absorption,  and  L (  e  )  is  the 
interaction  length  of  the  X-ray  and  prODellant 
sample.  The  distribution  of  X-ray  intensity 
continuous,  but  with  two  cusps  at  the  tanger. 
points  of  the  internal  perforation.  Even  thoug- 
the  existence  of  cusps  is  helpful  in  determining 
the  inner  diameter,  the  continuity  of  intensity 
profile  makes  it  difficult  to  use  direct  scaling 
and  standard-edge  detection  in  determining  the 
precise  internal  diameter  of  a  hollow  propellant 
sample.  Two-dimensional  slab  samples  woulc 
improve  the  situation  since  there  are  drastl : 
intensity  dlsoontinutlties  at  the  edges  between 
gas  and  solid  propellant  surfaces. 

Figure  p  snows  an  image  of  a  cylindrica. 
propellant  sample  and  an  X-ray  intensity  profile 
of  one  scan  line  perpendi cular  to  the  axi9 
obtained  in  an  actual  test  firing.  In  this 
picture,  the  valleys  of  intensity  are  not  aligned 
with  the  inner  diameter  due  to  the  scale 
difference.  Comparing  Fig.  6  with  Fig.  5,  it  is 
evident  that  the  image  obtained  in  the  test  firing 
is  degraded  severely  and  the  cusps  are  not  as 
distinct  as  those  in  Fig.  5.  The  physical 
mechanisms  responsible  for  this  degradation  are 
caused  by  l)  the  finite  focal-spot  size  of  the 
X-ray  tube  head,  which  generates  a  grainy  or 
unclear  image,  2)  Compton  scattering16  or  the  test 
chamber  and  propellant  sample,  which  becomes  more 
dominant  than  photoelectric  absorption  when  X-ray 
radiation  energy  becomes  higher,  3)  the  nonuniform 
intensity  distribution  of  X-ray  ’-adlatlon  over  the 
input  screen  of  the  Image  lntenslfler,  and  k)  the 
degradation  caused  by  the  image  intensifler, 
camera,  and  recording  system. 


Hwimtol  rMltln 

Fig.  6  X-ray  image  and  intensity  profile  of  a 
stick  propellant  sample 


To  accurately  deduce  the  Internal  dalmeter  of 
the  propellant  sample  by  Image  restoration,  the 
degradation  factors  mentioned  above  must  be 
properly  treated.  To  bypass  a  number  of  the 
difficulties  involved  in  image  restoration,  the 
following  procedures  are  adopted  in  determining 
the  instantaneous  location  of  the  internal  burning 
surface  of  the  propellant  sample. 

i  )  Generation  of  calibration  image:  A  set  of 
cylindrical  propellant  samples  were  machined  to 
specified  internal  diameters  for  calibration 
purposes.  The  calibration  images  were  produced  by 
placing  the  machined  propellant  samples  inside  the 
erosive-burning  test  rig  and  exposing  them  to 
X-ray  radiation  at  the  same  intensity  level  as 
test-firing  conditions.  The  setting  on  the  data 
acquisition  system  for  calibration  is  also 
identical  to  that  used  in  the  actual  test  firing. 
To  reduce  the  noise  level  in  the  calibration 
signal,  numerous  calibration  images  were  obtained 
and  averaged  for  a  single  specimen. 
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2)  Comparison  of  images  obtained  in  actual  test 
firing  with  calibration  Images:  The  intensity 
distribution  of  the  averaged  calibration  signals 
was  compared  with  the  intensity  distribution  of 
test-firing  images.  In  comparison,  a  normalized 
cross-correlation  given  by 


I  X(n)Y(n) 


is  used  where  X(n)  and  Y(n)  are  intensities  of 
X-ray  radiation  of  calibration  image  and  actual 
test  firing  image,  respectively. 

3)  Determination  of  Internal  diameter:  The 
instantaneous  internal  diameter  of  the  test 
propellant  grain  at  a  given  location  of  the  test 
propellant  grain  is  determined  from  the  best 
correlation  with  one  of  the  27  calibration 
propellant  samples.  Using  this  procedure,  the 
accuracy  of  the  diameter  measurement  is  in  the 
order  of  0.9  mm. 

Discussion  of  Results 

A  typical  set  of  pressur e- t i me  traces 
obtained  from  the  erosive-burning  test  firing  is 
shown  in  Fig.  7.  All  gages  show  pressure  spikes 
at  initial  time  less  than  0.1  s.  The  pressure 
spike  is  caused  by  the  strong  flow  of  product 
gases  generated  from  the  driving  motor.  After  the 
initial  spike,  the  pressures  at  all  gage  locations 
remain  relatively  constant  with  slight  decay. 


Fig.  7  A  typical  set  of  pressure-time  trac 
(Test  No.  ERB-29) 

X-ray  images  of  the  downstream  portion  of  the 
test  propellant  grain  were  recorded  on  a  video 
cassette  recorder.  (Due  to  the  length  limitation 
of  the  input  fluor  of  the  image  lntensifier,  only 
the  downstream  portion  is  recorded.)  A  typical 
X-ray  image  is  shown  in  Fig.  8.  I t  is  evident 
that  the  inner  diameter  of  the  propellant  grain 
can  be  identified.  The  increase  of  internal 
diameter  with  respec t  to  time  can  also  be  seen  on 
the  video  screen. 


Fig.  8  A  typical  X-ray  image  of  stick  propellant 
during  test  firing  ERB-29 


During  the  entire  test  period,  the  pressure  at  the 
driving  motor  is  higher  than  that  at  the  inlet 
nozzle  or  downstream  cavity.  This  implies  that 
the  flow  is  unidirectional  from  the  head  end  to 
the  exit  nozzle.  From  the  magnitude  of  the 
oressure  at  the  G3  location  one  can  be  sure  that 
except  at  the  end  of  the  test  firing,  the  exit 
nozzle  is  choked  most  of  the  time. 


After  comparison  with  calibration  samples  and 
carrying  out  data  reduction  procedures,  the  time 
variations  of  the  inner  diameter  distributions 
were  determined  (see  Fig.  9).  Using  the  data 
presented  in  Fig.  9,  the  burning  rates  at  the 
upstream  (x*xu*6.03  cm)  and  downstream  (x-xd-15.58 
cm)  locations  were  deduced  and  plotted  in  Fig.  10. 
The  strand  burning  rate  of  the  NOSOL-363 
propellant  was  also  calculated  for  evaluation  of 
eroslvlty  (rb/rb°)  at  the  two  locations  described 
above.  This  information  is  plotted  together  with 
the  deduced  burning  rate  on  the  same  figure. 
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•'ig.  9  Time  variation  of  the  Inner  diameter 
distribution  deduced  from  recorded  X-ray  Images. 


Fig.  10  Measured  time  variation  of  eroslsve- 
burnlng  rates  and  erosivlties  at  upstream  and 
downstream  locations. 


Some  interesting  observations  made  from  Figs. 

9  and  10  are  summarized  below: 

1)  In  the  early  time  period  C t < 1 3 ) ,  the  con¬ 
tour  of  the  inner  diameter  exhibits  diver¬ 
gence  in  the  downstream  direction  due  to 
the  strong  erosive  burning  effect  exhibit¬ 
ed  there.  The  flow  Mach  number  near  the 
downstream  end,  xp,  is  nearly  sonic,  and 
hence  has  very  high  erosivity  in  the 
order  of  3.2.  The  Mach  number  near  the 
xu  location  is  also  quite  high  (-0.6), 
but  the  erosivity  i3  about  half  of  that 
at  the  xD  location.  During  this  time 
period,  the  pressure  gradient  (-  3p/3x) 
decreases  drastically  with  respect  to 
time,  and  the  perforation  area  increases 
rapidly.  Both  effects  cau3e  the  erosiv¬ 
ity  at  the  Xp  to  decrease. 

2)  Between  tj  and  tj,  the  burning  rate  at 
the  upstream  location  catches  up  to  and 
exceeds  that  at  the  xp  location.  This  is 
caused  by  the  fact  that  as  the 
perforation  area  opens,  the  flow  Mach 
number  decreases  and  the  cross-flow 
effect  becomes  less  pronounced.  Due  to 
the  decrease  in  pressure,  however,  the 
erosivity  may  even  exhibit  a  peak  region 
before  rapidly  decaying.  As  a  result  of 
the  higher  burning  rate  at  xu,  the 
contour  of  the  internal  perforation 
becomes  uniform. 

3)  Between  tj  and  t 7,  the  burning  rate  and 
erosivity  both  decrease  with  respect  to 
time  at  all  locations.  During  this 
period,  the  propellant  charge  In  the 
driving  motor  is  largely  consumed. 
Because  the  mass  flow  rate  In  the 
crossflow  is  produced  mainly  from  the 
propellant  sample,  the  mass  flow  rate  in 
the  downstream  location  becomes  higher 
than  that  of  upstream.  Thus,  the  rate  of 
decay  of  the  burning  rate  at  xp  is  slower 
than  that  at  xu.  As  a  result  of  this 
variation  in  burning  rate,  the  contour  of 
the  internal  perforation  again  becomes 
divergent. 

A)  The  diameter  distribution  of  the 
recovered  propellant  grain  between  xu  and 
xp  is  shown  in  Fig.  9.  This  information 
was  obtained  from  the  X-ray  image 
presented  i  r.  Fig.  11.  The  contour  of 
the  Port  area  was  quite  interesting.  The 
minimum  web  thickness  occurred  at  x  - 
6.05  cm  and  a  portion  of  the  grain  even 
exhibited  a  burn-through .  The  X-ray 
image  shows  a  divergent  section  in  the 
head  end  of  the  propellant  grain;  this  is 
due  to  the  design  of  the  inlet  graphite 
nozzle  which  protrudes  into  the 
propellant  grain.  Near  the  end  of  test 
firing,  as  other  portions  of  the 
propellant  surface  are  extinguished  by 
the  low  pressure  condition,  the  weak 
gaseous  Jet  from  the  inlet  nozzle  can 
still  expand  and  impinge  on  the  inner 
surface  of  the  propellant  grain. 
Therefore,  the  burning  in  this  portion 
during  the  extinction  period  produced  a 
convergent-divergent  contour. 
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Fig.  11  X-ray  Image  of  recovered  propellant  grair  from  teat  firing  ERB-24 


Summary  and  Conclualona 

1)  A  comprehensive  theoretical  model  has  been 
formulated  to  simulate  er os  1 v e- bur n i n g 
processes  of  stick  propellants. 

2)  The  feasibility  of  using  real-time  X-ray 
radiography  to  measure  the  Instantaneous 
propellant  burning  surface  has  been 
demonstrated. 


Turbulent  Reacting  Flows,  edited  by  Libby,  P. 
A.  and  Williams,  F.  A.,  Sprlnger-Verlag , 
Berlln-Heldelberg-New  York,  1980. 

Janicka,  J.  and  Kollman,  W.,  Seventeenth 
Symposium  (International)  on  Combustion , 
pp.  42 1-430,  1979. 

Bray,  K.  N.  C.  and  Moss,  J.  B. ,  Combustion 
and  Flame.  Vol.  30,  pp.  125-131.  197^ 


3)  Based  upon  the  recorded  X-ray  images,  the  11. 
Instantaneous  burning  rate  of  NOSOL-363  stick 
propellant  has  been  determined  under  test 
motor  operating  conditions.  Results  show  the  12. 
very  strong  influence  of  cross-flo*  velocity 
on  propellant  burning  rate.  (Eroslvlty  can  be 
as  high  as  3.2.)  This  implies  that  in  order  13. 
to  achieve  accurate  and  realistic  predictions 
of  stick-propellant  combustion  performance  in 
gun  propulsion  systems,  the  erosive-burning 
effect  must  be  properly  Incorporated.  14. 
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